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Abstract 


The Swift/XRT detected the X-ray afterglow of long burst GRB 2201174, which began to rebrighten 300 s after 
triggering and followed a single power-law decay segment after thousands of seconds of the orbital observation 
gap. This segment is different from the shallow decay segment (plateau) and flare, and may belong to a giant X-ray 
bump. We investigated this segment by the fall-back accretion model and found that the model can interpret this 
segment with reasonable parameter values. Within this physical model scenario, the fall-back accretion rate reaches 
a peak value ~1.70 x 10 ^M, s^! around 300s in the central engine frame, which is compatible with the late 
mass supply rate of some low-metallicity massive progenitor stars. The initial black hole (BH) spin is 
ao = 0.647024 and implies that this re-brightening signature requires a larger BH spin. The total accretion mass 


during the fall-back process is Mace = (3.09 + 0.02) x 107M. 


The jet energy from the fall-back accretion is 


(9.77 + 0.65) x 10?? erg, with a ratio of 0.066 to the isotropic-equivalent radiation energies of the GRB prompt 
phase in the 1-10* keV band. The fall-back radius Tp corresponding to the peak time of fall-back t, is 
(3.16 + 0.05) x 10!? cm, which is consistent with the typical radius of Wolf-Rayet stars. In summary, our results 
provide additional support for the origin of the long burst from the core collapse of Wolf—Rayet stars, and its late 
central engine activity is likely due to the fall-back accretion process. 
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1. Introduction 


Gamma ray bursts (GRBs) are the brightest electromagnetic 
event that have occurred in the universe since the Big Bang, 
consisting of two emission phases: the prompt emission(with 
an initial prompt soft y-ray emission) and afterglow emission 
(with long term broadband emission). The prompt emission is 
generally considered to be related to the internal dissipation of 
the jet, such as internal shock dissipation or magnetic 
dissipation. Afterglow emission is usually considered to be 
from external shocks (especially forward shock) generated by 
the interaction between jets and interstellar medium (see 
Zhang 2018 for a review). In general, the end of the prompt 
emission phase means that the cease of GRB’s central engine. 
However, observations from Swift satellite indicate that many 
GRB’s central engines have an extended activity time, mainly 
manifested as shallow decay segments (plateaus) (Liang et al. 
2007; Troja et al. 2007; Tang et al. 2019; Zhao et al. 2019), 
flares (Burrows et al. 2005; Zhang et al. 2006; Margutti et al. 
2011), and giant bumps (Kumar et al. 2008a; Wu et al. 2013; 
Gao et al. 2016; Chen et al. 2017; Zhao et al. 2020, 2021) in the 
X-ray light curves following the prompt emission. 

The Swift/XRT detected the X-ray afterglow of long burst 
GRB 220117A, which began to rebrighten 300s after 
triggering and suggest that this GRB has an extended central 
engine activity time. In addition, this signature is different from 
X-ray shallow decay segment (plateau) and flare, and may 


belong to a giant X-ray bump. Due to the absence of an 
“internal plateau” feature in the X-ray afterglow of GRB 
220117A, its central engine may be a hyperaccreting black hole 
(BH) system. In this physical scenario, the internal dissipation 
of fall-back accretion energy can interpret the giant X-ray 
bump. If the fall-back accretion rate or the duration of fall-back 
accretion process is large enough, a giant X-ray bump that 
rapidly rises and decays with a form of f° /3 is expected. So 
far, the giant X-ray bump has been found in the X-ray 
afterglow of many GRBs and could be well interpreted within 
the fall-back accretion model (Wu et al. 2013; Gao et al. 2016; 
Chen et al. 2017; Zhao et al. 2021). 

In this paper, we study the X-ray re-brightening signature of 
GRB 220117A within the fall-back accretion model. In 
Section 2, we describe the observations of GRB 220117A 
and analyze the re-brightening signature. The fall-back 
accretion model is described in Section 3. In Section 4, we 
apply the fall-back accretion model to the re-brightening 
signature. The conclusion and implications of our results are 
discussed in Section 5. Throughout the paper, the convention 
Q = 10"Q, is adopted in c.g.s. units. 


2. GRB 220117A Observations 


The BAT triggered and located GRB 220117A at 23:58:21 
UT on 2022 January 17. T99(15—350 keV) is 49.81 + 2.37 s. 
The time-averaged spectrum from T+ 14.79s to T+ 66.14s is 
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Figure 1. The XRT light curve of GTB 220117A. The solid red lines are the best fits with a smooth broken power law for the X-ray re-brightening phase and its 


follow-up decay phase. 


best fitted by a single power law (SPL) function. The power- 
law index of the time-averaged spectrum is [,=1.8+ 0.18. 
The fluence in the 15-150keV band is $,=1.620.2 x 
1076 erg cm "(Palmer et al. 2022). XRT observations started 
at 151.9 s after the trigger. The XRT light curve is shown in 
Figure 1. The UVOT starts collecting data 161s after the 
trigger. No source was detected by the UVOT at the X-ray 
afterglow position (Melandri et al. 2022). Palmerio et al. (2022) 
observed the afterglow of GRB 2201174 using the ESO VLT 
UT3. From the feature of Lya trough and the Sill 1260 AA, 
the redshift was measured as z — 4.961. 

For the re-brightening segment and follow-up segment of 
GRB 220117A, we adopted a smooth broken power law (BPL) 
function to fit it, 


t way t Wa -1/s 
ep o 
fp Ip 


where o, and o» represent the decay scope before and after the 
break, respectively. w describes the sharpness of break, and 
here we adopt w — 3 as suggested by Liang et al. (2007). The 
fiting light curve and the best fitting parameters are shown in 
Figure 1. In 2019, Zhao et al. (2019) found that the decay scope 
of a shallow decay segment accords with normal distribution 
ay — 0.35 £0.35 based on 13 yr of Swift/XRT observation 
data. It can be inferred that X-ray re-brightening segment may 
not be a shallow decay segment. 


In addition, we also tested whether the X-ray re-brightening 
segment and follow-up segment belong to the flare. Yi et al. 
(2016) analyzed the GRBs with significant flares observed by 
Swift/XRT from 2005 April to 2015 March and obtained an 
empirical relationship: 


logi Taur = (—0.35 + 0.11) + (1.12 + 0.04) x logyo (Ts) 
Q) 


where Tpeak is the peak time of the fare. Taur = Tena — Tstart 18 
the duration of the flare, where 7.444 and Tena represent the start 
time and end time of the flare, respectively. Tpeak and Taur of 
each flare could be easily obtained by fitting the light curve 
with a smooth BPL function. We find that the duration of X-ray 
re-brightening segment and its follow-up decay segment 
deviates from the duration obtained by the empirical relation- 
ship 2c. Compared to the typical flare, the giant X-ray bump 
has a relatively longer duration. Therefore, the X-ray re- 
brightening segment and follow-up segment may belong to a 
giant X-ray bump. 


3. Model Description 


In this paper, we intend to use the fall-back accretion model 
to interpret the X-ray re-brightening segment of GRB 
220117A. The physical scenario of the fall-back accretion 
model is described as follows: the progenitor stars of long 
GRBs may be massive stars (Woosley 1993; Paczyński 1998; 
MacFadyen & Woosley 1999; Woosley & Bloom 2006), and 


Research in Astronomy and Astrophysics, 23:115014 (7pp), 2023 November Zhao 


generally have a core-envelope structure. At the end of the 
massive stars' evolution, photodisintegration and electron 
capture will trigger core collapse and form a super accretion 
BH system. The relative jet is launched by the super accretion 
BH system through extracting the gravitational energy from the 
accreted material or the spin energy of the BH. When the 
relativistic jet successfully penetrates the envelope, it powers 
initial prompt emission and broadband afterglow emission of 
the GRB. During the process of jet penetration through the 
progenitor envelope, a part of jet energy is transferred to the 
envelope, which might help the supernova to explode. The 
bounding shock responsible for the associated supernova 
transfers kinetic energy to the envelope materials and a part 
of envelope materials would be ejected. The rest of envelope 
materials could fall back into around of BH and form an 
accretion disk (Kumar et al. 2008a, 2008b). The accretion disk 
may power a new relativistic jet through the neutrino- 
annihilations mechanism (Popham et al. 1999; Narayan et al. 
2001; Di Matteo et al. 2002; Janiuk et al. 2004; Gu et al. 2006; 
Chen & Beloborodov 2007; Liu et al. 2007, 2015; Lei et al. 
2009, 2017; Xie et al. 2016) or Blandford-Znajek (BZ) 
mechanism (Blandford & Znajek 1977; Lee et al. 2000; 
Li 2000; Lei et al. 2005, 2013). Compared with the neutrino- 
annihilations mechanism, the jet powered by the BZ mech- 
anism will be cleaner and more powerful (Lei et al. 2017; Xie 
et al. 2017; Lloyd-Ronning et al. 2018). Therefore, the BZ 
mechanism is more likely to interpret the late central engine 
activity of long GRBs. A part of BZ energy would eventually 
be injected into the afterglow blast wave. If the injected energy 
is comparable to or even larger than the blast wave kinetic 
energy, it will generate some detectable signature such as a 
plateau in the X-ray afterglow (Zhao et al. 2020). In addition, 
the rest energy would undergo the internal dissipation process, 
which may generate the observed giant X-ray bump. 

The evolution of the fall-back rate of progenitor envelope 
material can be described by a smooth BPL function as 
(Chevalier 1989; MacFadyen et al. 2001; Zhang et al. 2008; 
Dai & Liu 2012) 


(ei) I cmi ie 
Ma = M,| -| —* 4p E 
T [A | cer | Si 


where fg is the starting time of the fall-back accretion in the 
central engine frame, M, is the peak of fall-back rate, and 1, is 
the peak time of fall-back rate in the central engine frame. 
Progenitor envelope material falls back in the late time and 
forms an accretion disk, which will be viscous accretion by the 
BH. The accretion rate can be obtained as (Kumar et al. 20082) 


woes (4) 


Tyis 


where Tyis ~ 1/aQ, is viscous timescale of accretion disk. a is 
a standard dimensionless viscosity parameter with values of 


~0.01-0.1, and €), is the Kepler angular velocity of the 
accretion disk. Due to accretion disk mass M, will increase 
with the fall-back from the envelope material and decrease with 
accretion. Therefore, (Kumar et al. 20082; Lei et al. 2017) 
Ma = My — M. (5) 


By combining Equations (4) and (5), the accretion rate can be 
obtained as 


M= 


t 
1 etis f es Madt', (6) 
Tyis to 

If the viscous timescale is larger than the accretion time, the BH 
will undergo a slow accretion process. The accretion rate of BH 
would be flat when ¢>7,;, and then begin to decrease with 
time, see Figure 7 in Lei et al. (2017). On the other hand, the 
BH will undergo rapid accretion process if Tyi «t. The 
evolution of the accretion rate will follow the fall-back rate, i.e., 
M = Mp. In this paper, we assume that the accretion process of 
BH undergo a rapid accretion process. 

The fall-back accretion can continuously extract the 
rotational energy of the BH through the BZ mechanism, and 
powers a Poynting dominated jet. A part of jet energy would 
undergo the internal dissipation process, which may generate 
the observed giant X-ray bump. To connect the observed X-ray 
luminosity Ly and BZ power Pgz, it is necessary to introduce 
efficiency factor ņy and jet beaming factor f, 


NxPpz = f, Lx- (7) 


The power of the BZ mechanism to extract the rotational 
energy of BH can be written as (Lee et al. 2000; Li 2000; Wang 
et al. 2002; Lei et al. 2005, 2013, 2017; McKinney 2005; Lei & 
Zhang 2011; Chen et al. 2017; Liu et al. 2017; Lloyd-Ronning 
et al. 2018) 


Pgz = 1.7 x 10°a2m? B25 f(a.) erg s^, (8) 


with 


2 
q=Vl-a, (9) 


jajsi- (=). 


where a. is the dimensionless BH spin parameter, m. = M./Mo 
is the dimensionless BH mass, and B.;5 is strength of the 
magnetic field near the BH horizon in units of 10? G. From the 
Equations (8) and (9), it can be seen that the BZ power mainly 
depends on parameters of B., m. and a.. 

In general, the strength of the magnetic field near the BH 
horizon can be estimated by balancing the magnetic pressure 
on the BH horizon and ram pressure of the accretion flow at the 
inner edge of the accretion disk (Moderski et al. 1997) 


B2 " Mc 
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Figure 2. The fitting result of GRB 220117A’s X-ray re-brightening segment by adopted the fall-back accretion model. 


where n = (1 + y1 — a23)GM./ c? is radius of BH horizon. 

The BZ process extracts rotational energy and angular 
momentum from the BH, while the accretion process brings the 
energy and angular momentum of the accretion disk into the 
BH. According to the conservation of energy and angular 
momentum, the evolution of BH under the two processes can 
be written as (Wang et al. 2002) 


2 
en = Mc Ew z Pgz, (11) 
dt 
eh = LinsM E Tez. (12) 
dt 


From above two equations, it can be derived that 


da.  (MLm — Tsz)c — 2a.(Mc!?Es, — Pgz) 
dt GM? M.c? 


» (13) 


where Tz is the torque applied to the BH by the BZ process, 
which can be written as (Bardeen et al. 1972) 


Tez = 3.36 x 10%a?q m2 B2 sF (a)g cms. (14) 
Here Ems and Lms are the specific energy and angular 
momentum at the radius of the innermost inner edge of 
accretion disk Risco, respectively, defined as (Novikov & 


Thorne 1973): 


Ems = (44 Risco = 3a.) / (43 Risco), (15) 


Lms = (GM./ c) (2 (3X Risco xm 2a)/(4 3R ico); (16) 


with 


TEB ta- IG -a6 +a + 227], 


za=1 +4 0 - aP + a)! 4 (1 — 4), 
z = (Ba? + 2)”. 


Risco(M., a.) ms 


(17) 


4. Model Application to the X-Ray Re-brightening 
Signature of GRB 220117A 


In this section, we apply the fall-back accretion model that 
was introduced in Section 3 to the X-ray re-brightening segment 
and follow-up segment of GRB 220117A. We adopt the start 
time of re-brightening segment faa = fstart,obs/ (1 +z) and the 
end time of follow-up segment fend = fendobs/(1 + z) as the start 
time and the end time of the fall-back accretion in the central 
engine frame. According to analysis, the initial mass of BH 
hardly affects the BZ power. In this section, we adopt 
M.9— 3M... In addition, we take f; — 0.01 and 5y — 0.01 in 
our calculation. Finally, we take the dimensionless peak fall- 
back rate m,, the sharpness of the peak s, the peak time of fall- 
back rate t, and initial BH spin ao as our free parameters. In 
order to obtain the best fitting values, a Markov Chain Monte 
Carlo (MCMC) method is adopted. In our MCMC fitting, the 
emcee code is used (Foreman-Mackey et al. 2013). In the 
code, we set the boundaries of four free parameters to 
logio) = [—15, 0], s=[0, 10], tp = [tse fena] and ao = [0, 
1], respectively. 

The fitting result of GRB 220117A's X-ray re-brightening 
segment and follow-up segment is shown in Figure 2. In 
addition, Figure 3 shows the corner plot of the free parameters 
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Figure 3. The corner plot of the free parameters posterior probability distribution for the fitting result. 


posterior probability distribution. We adopt the median of the 
free parameter distribution and lo error as the fitting results, 
with the dimensionless peak fall-back rate is logo) = 


ATT Os, the sharpness of the peak is s = 0.60* 010. the peak 


time of fall-back rate is log,)(t,) = 2.51 + 0.02 and the initial 
BH spin is ag = 0.647024. From the fitting results, we 
calculated the total mass during the fall-back process is 
Mace = (3.09 + 0.02) x 10 ?M., the strength of the magnetic 
field near the BH horizon at the ¢, is Bj,,5 = 0.089 + 0.0021, 
the fall-back radius corresponding to t, is r, = (3.16 + 0.05) x 
10'°cm and the jet energy from the fall-back accretion is 
(9.77 £ 0.65) x 10°erg, with a ratio of 0.066 to the isotropic- 
equivalent radiation energies of the GRB prompt phase in the 
1-10* keV band. 

In conclusion, the fall-back accretion model can interpret the 
X-ray re-brightening segment and follow-up segment with 
reasonable parameter values. This signature may require a 
larger black hole spin (the peak of its posterior probability 
distribution is ay = 0.64^026). When the fall-back rate reaches 
its peak, the corresponding fall-back radius is consistent with 
the typical radius of a Wolf Rayet star. This result provides 
further support for the origin of the long GRB from the core 
collapse of the Wolf-Rayet stars, and its late central engine 
activity is likely due to the fall-back accretion process. 

The fallback accretion rate reaches a peak value 
~1.70 x 10°-°M, s^! around 300 s in the central engine frame. 
The bounding shock responsible for the associated supernova 


transfers kinetic energy to the envelope materials and most of 
the envelope materials would be ejected. Therefore, it is very 
necessary to check whether the progenitor of GRB 220117A 
can provide enough envelope material. The fall-back rate of 
envelope material can be estimated by the data of pre-SN 
models, i.e., (Suwa & Ioka 2011; Woosley & Heger 2012; 
Matsumoto et al. 2015; Liu et al. 2018) 


My = CES } (18) 
tee Vp — p 


where p is the average density within radius r, t; is the freefall 
timescale of envelope material with radius r. It can be obtained 


as 
uM Ea (19) 
" 2 \ 2GM, 


where M, is the mass within the radius r. The value can be 
obtained by the following equation 


M, = M.o + f 4nr?pdr, (20) 
ro 


where p is the density of envelope material at radius r. We set 
the time and radial coordinate for which the enclosed mass 
reach initial black hole mass as rọ and 0, respectively. By adopt 
the progenitor density profile with different metallicities and 
masses from Liu et al. (2018), we calculate evolution of the 
mass supply rate for theirs and are shown in Figure 4. We find 
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Figure 4. The evolution of the mass supply rate of progenitor with different metallicities and masses. The dashed line represents fall-back rate evolution of GRB 
220117A's progenitor envelope material in our fitting. Zo is the metallicities of the Sun. 


that low-metallicity massive progenitor stars are compatible 
with our fitting results. 


5. Discussion and Conclusions 


The Swift/XRT detected the X-ray afterglow of long burst 
GRB 220117A, which began to rebrighten 300 s after 
triggering and followed a single power-law decay segment 
after thousands of seconds of the orbital observation gap. The 
re-brightening segment is different from the shallow decay 
segment (plateau) and flare, and may belong to a giant X-ray 
bump. We investigated this segment by the fall-back accretion 
energy internal dissipation model. We found that the model can 
interpret this segment with reasonable parameter values. Within 
this physical model scenario, the fallback accretion rate reaches 
a peak value ~1.70 x 10 ^M. s^! around 300s in the central 
engine frame, which is compatible with the late mass supply 
rate of some low-metallicity massive progenitor stars. The 
initial BH spin is ao = 0.64*02« and implies that this re- 
brightening signature requires a larger BH spin. The total 
accretion mass Mc during the fall-back process is 
Mace = 3.09 x 10 7M... The jet energy from the fall-back 
accretion is (9.77 + 0.65) x 10° "erg, with a ratio of 0.066 to the 
isotropic-equivalent radiation energies of GRB prompt phase in 
the 1-10* keV band. The fall-back radius rj corresponding to t, 
is (3.16 + 0.05) x 10!°cm, which is consistent with the typical 
radius of Wolf-Rayet stars. In summary, our results provide 
additional support for the origin of the long burst from the core 
collapse of Wolf-Rayet stars, and its late central engine activity 
is likely due to the fall-back accretion process. 

In this paper, we adopted a simple fall-back accretion rate 
evolution model, and do not consider the angular momentum 


distribution of the progenitor star, so this calculation is 
approximately valid for a slowly rotating progenitor. In fact, 
the angular velocities €? of stars with different radii are different 
and the fall-back radius satisfies x9? (Kumar et al. 20082). 
Therefore, the angular momentum distribution of the progenitor 
star has a great influence on the fall-back accretion rate. In 
addition, the metallicities and masses of stars will have a 
certain impact on the fall-back accretion rate of the envelope 
material. In the future, we will study the fall-back accretion 
model with different angular momentum distribution, mass 
distribution and metallicity. 

Super Eddington accretion makes material on the accretion 
disk subject to an outward radiation pressure that is greater than 
gravity, Therefore, an outflow (i.e., disk wind) driven by 
radiation pressure is launched on the surface of the accretion 
disk, and taking away a part of the fall-back materials. 
However, the disk wind is ignored in the fall-back accretion 
model adopted in this paper. It will cause the accretion rate of 
BH to decrease. Because the outflow of the accretion disk is not 
well understood, the power law function model is generally 
adopted to describe the accretion rate at different accretion disk 
radii. It can be seen that the influence from the outflow of the 
accretion disk to the accretion rate is highly dependent on the 
power law index. In addition, the existence of the accretion 
disk outflow will also be important for understanding the 
baryon load of a GRB jet (Lei et al. 2013, 2017) and 56Ni 
synthesis for associated supernovae (Song & Liu 2019). We 
hope that the future general-relativistic magnetohydrodynamic 
(GRMHD) simulations can help us better understand the 
accretion disk outflow. 
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